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(54) A LINEAR INTERPOLATING METHOD FOR COLOR 
SIGNALS IN A MEMORY 

(71) We, DAINIPPON SCREEN SEIZO KABUSHIKJ KAISHA, a 
Japanese Corporation, of 1-1 Tenjin-kitamashi, Teranouchi-agaru 4-chome, Horikawa- 
dori, Klamigyo-ku, Kyoti-shi, Japan, do hereby declare the invention, for which we 
pray that a patent may be granted to us, and the method by which it is to be per- 
5 formed, to be particularly described in and by the following statement:-— 5 

This invention relates to a linear interpolating method for signals in a memory 
which is used for color correction of picture signals in a reproducing machine such 
as a color scanner, a color facsimile producer, or the like, in which color separation 
picture images are produced by photo-electric scanning. 
10 in conventional color photographic plate making, color correction is often made jq 

by. photographic masking. However, this method has many defects, for example: 
limitations of color correction ability, necessity for many skilled engineers, unreliable 
results of the color separation, irregular quality of finish, complexity, and the like. 

In order to overcome these defects, a color correction masking method by an 
15 electronic color separation machine such as a color scanner has been developed and is 15 
nowadays more popular. Most of the color scanners now used employ an analog 
computer system for the color correction calculations so as to increase the calculation 
speed. 

This method, however, has also defects such as the difficulty of the introduction 

20 of many kinds of calculations because of the restriction of calculation ability, inevitable 20 
effects of temperature drift and noise, multiplicity of operational amplifiers and so 
forth as electric elements, inconvenience of operation due. to many adjustments of 
potentiometers and switches, and high manufacturing cost. 

If the analog computer system is simply replaced with a digital computer system, 

25 which has advantages such as a wide correction variable range and convenience of 25 
operation, the calculation speed for the color correction decreases very much, and the 
processing ability is reduced. Accordingly, this system is not practicable. 

Recently, a direct scanner has been developed for plate making in printing, which 
performs color separation, color correction, conversion of scale of the reproduced 

30 image, and halftone processing at the same time so as to meet the requirement for 30 
high quality printing and rapid operation. In this case, however, there is. the defect 
that supplementary masking or hand retouching after the color separation cannot be 
applied, as opposed to conventional color scanning which includes color separation, 
color correction, conversion of scale of the reproduction image, and halftone processing. 

35 In general, an original color picture is scanned by a color scanner to obtain three 35 

(red, green, and blue) color separation signals. These three color separation signals 
are sent to a color operation circuit, thereby finally obtaining recording signals for 
density of printing inks, such as cyan, magenta, yellow, and black. 

In. order to provide the most accurate possible color reproduction, a combination 

40 of the amounts of cyan, magenta, and yellow inks (the black ink, and so forth, are 40 

omitted for the sake of brevity of explanation) is necessarily determined corresponding 
to a combination of red, green, and blue color separation signals. 

Consequently, for the purpose of color correction by selecting the combination 
of cyan, magenta, and yellow values corresponding to the combination of red, green, 

45 and blue values/ the color-corrected combinations of cyan, magenta, and yellow values 45 

corresponding to each combination of red, green, and blue values are stored in a 
memory in advance, and then the color-corrected combination of cyan, magenta, and 
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yellow values is read out by addressing the memory by the combination of red, green, 
and blue values corresponding thereto. 

If each red, green, and blue range is divided into, for example, two hundred tone 
steps, altogether 200 s =8,000,000 combinations of cyan, magenta, and yellow values 
5 must be stored in the memory, which requires that the memory have a large capacity. 

This means high cost, and thus is not practicable. .7 

Therefore, in order to reduce the storage capacity required for the memory, each 
color range of red, green, and blue is divided into, for example, sixteen tone steps, 
and then 16* = 4096 combinations of cyan, magenta, and yellow values are required. 
10 Thus the storage capacity requirement for the memory is reduced to a manageable jq 
level. On the other hand, the tone steps become too rough, and the lack of output 
consistency becomes conspicuous, so that printing quality suffers. Therefore, in this 
case, it is necessary to interpolate intermediate values properly between each two tone 
steps. 

15 The present invention seeks to provide a linear interpolating method for such 15 

ssignals in a memory free from the above mentioned defects, which enables the ' 
memory to calculate quickly by using a simple formula interpolation values, without 
large discontinuities of the slope of the interpolated values between one interpolation 
unit and the next, 

20 According to the present invention there is provided a linear interpolating 20 

method for color signals in a memory of a picture reproducing machine,, comprising 
storing appropriate values of color picture output signals corresponding to certain 
stepped values of color input signals in the memory addressed in a three-dimensional 
fashion, and interpolating values of color output signals at points which are between 

25 said values by dividing up the cubic interpolation unit of the memory which is 25 

constituted by a single step of each of the color input signals into a plurality of 
tetrahedra whose vertices are either vertices of the cubic unit, centers of its faces, 
or its center point, calculating the color output signal at each vertex of these tetra- 
hedra which is a center of a face of the cubic unit, if any, by averaging the values 

30 of the color output signal at the four vertices which are corners of said face, and 30 
at the center point of the cubic unit by averaging the values of the color ouput signal 
at all eight of the vertices of the cubic unit, determining which of these tetrahedra 
includes the interpolation point at which the value of the color output signal is to be 
interpolated, and deriving the interpolated value at the interpolation point as a 

35 weighted sum of the values at the four vertices of the determined tetrahedron, the 35 
value at each vertex being given a weight corresponding to the ratio of the volume 
of a second tetrahedron whose vertices are the interpolation point and the other three 
vertices of the determined tetrahedron to the volume of the determined tetrahedron. 
In one embodiment in accordance with the invention the cubic unit is divided 

40 into twenty-four tetrahedra, of each of which one vertex is the center of the cubic unit, 40 
two vertices are vertices of the cubic unit which are connected by an edge of the cubic 
unit, and the fourth vertex is the center of a square face of the cubic unit one edge 
of which face is the said edge of the cubic unit. 

In a second embodiment the cubic unit is divided into six tetrahedra by three 

45 planes which have a line in common, said line being a long diagonal of the cubic unit, 45 

• and each plane containing two edges, and four vertices of the cubic unit. 

In a still further embodiment the cubic unit is divided into five tetrahedra by 
four planes, each of which contains exactly three vertices of the cubic unit, said 
planes intersecting one another along lines which are diagonals of the faces of the 

50 cubic unit. ~ 50 

The present invention will now be described in detail with respect to the 
accompanying drawings, in which: 

Fig. 1 is a schematic view of a conventional interpolating method over a two- 
dimensional interpolation unit square; 

55 Figs. 2 and 3 are schematic views of a square interpolation unit region and a 55 

cubic interpolation unit region of the conventional two-dimensional and three- 
dimensional interpolating methods; 

Fig. 4 is a schematic view of a distribution of interpolation values of the conven- 
tional method for two-dimensional interpolation; 

60 Fig. 5 is a schematic view of an improved two-dimensional interpolting method; 60 

Fig. 6 is a schematic view of a method of three-dimensional interpolation over a 
tetrahedral region; 

Fig. 7 is a schematic view of a cubic interpolation unit dissected into twenty-four 
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tetrahedra according to one of the variations of the method of the present invention; 
Fig. 8 is a schematic view of one. of the tetrahedra of the dissection of Fig. 7; 
. Fig. 9 is a schematic view of another variation, of the method of the present 
invention, wherein the unit cube is dissected into six tetrahedra; 
5 . : : Fig. 10 is a schematic view of the tetrahedra obtained by the dissection of Fig. 9; 5 

and'. 

J Figs. 11 and 12 illustrate another method of dissecting the unit cube into five 
tetrahedra, which gives another variation of the method of the present invention. 

In order that the method may be better understood, some explanation of prior 
10 art methods of interpolation will now be given. 10 

Referring to Fig. 1, there is shown an example of interpolation of a function U 
of two variables, where the interval to be interpolated over is taken as unity. 

The value U(x,y) i.e. U(xj+Xr, yt+Vf) at a point P in an interpolation region 
ABCD will be found by a methematical interpolating method, in which x t and y^ are 
15 the integral parts of x and y and x, and y f are the decimal parts. 15 

For the interpolation it is necessary that the function at the vertices A, B, Q 
. and P should have known values U(x b yi), U(xi+l,yi), U(x t -f l,yj + l), and 
. yrx^Hbl). The interpolated value U(x,y) will be a function of x f , y f , U(xbVi), 
U(xj4-.i,y0, U(Xf-r*l,yi+l), and U(x b yi+1). Further, for consistency, the inter- 
20 polated value should be consistent with the known values of the original function at 20 
the corners of the unit region. . 

Ah interpolating method satisfying such a condition will be described. It is 
called linear interpolation because on the edges of the unit region it reduces to a 
simple linear interpolation function. 
25 In order to find the value U(x,y) at the point P in the interpolation unit square 25 

ABCD, first draw four perpendiculars from the point P to each side AB, BC, CD, 
and DA of the square. Designate the feet of these perpendiculars by Qj, Q 2 , Q 8 , and 
Q 4 respectively, as shown in Fig. 2, and add up the results obtained by multiplying 
each known value at the vertices A, B, C, and D by the area of each rectangle opposite 
30 to the vertex, thereby obtaining the following equation (I): 30 

Vfry^Vfr+XoK+yf^Vtw) (1-*) (!-*) 

+U(*+l,y0 x, (1-y,) , (I) 
+U(x„y l +1) (I-*) y, 
+U(*+l,yt+l).s* y, 

The interpolating method according to the formula (I) satisfies the above 
boundary conditions at the corners of the unit square and reduces to linear interpolation 
along the edges of the unit square, and thus is mathematically reasonable. Further, 
35 this method may be applied to the three-dimensional case. 35 

In Fig, 3 there is shown a unit cube interpolation unit having eight vertices 
with co-ordinates of 

Ota*)* (x,+l,yiZt), 
(x*,y,+ l,Zi), (xt^-f 1), 

40 (Si-Hltfi^- l^zi), («+layu%+l)j 40 

(xbyt+l^-hl), and (^-hl^+lA-f 1)* 

and including, a point P which co-ordinates (xi+x ti yi+y/, Zi+zi) at which the value 
of U is to be interpolated. The cube is divided up into eight rectangular paralello- 
pipedia by three planes which include the point P and are parallel to its faces. The 
45 value U(x,y,z) at the point P is found by adding up the values obtained by multiplying 45 
each known value at each of the vertices of the unit cube by the volume of each 
rectangular parallelopipedon which is positioned opposite to that vertex, thereby 
obtaining the following formula (II): 
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U(w>=U(*.+*t, y.+yr. 2 l +z / ) = U( 3 c i) y lA ) (1-*) (l^y,) (1-z,) 

+u(*+ijfo*) * (i— y.) 

+U(xur,+1*) (l-x t ) y, (l-.z.) 
+U(x,,y b 2 1 +1) (1-*) (l-y f ) z, 
+U(x„y,+U+1) (1-x,) y, z, 
+U(x,+l,y l ,z 1 +l) x, (l-y f ) z, 
+U(*+1*+U) * y, (1-z,) 
+U(x,+l,y<+U+l) x, y, z, 



(n) 



Again, this method produces consistent results at the vertices of the unit cube. 
Further, along the edges of the unit cube it reduces to simple linear interpolation, 
and on the faces of the unit cube it reduces to the method of quation (I). It is further 
5 clear that the value obtained in the center of each face of the unit cube is the mean ^ 
value of the known values at each vertex of that face, and the value obtained at the 
center of the unit cube is the mean value of the eight known values at the vertices 
of the cube. Accordingly, this method is seen to be mathematically reasonable. 

However, this method has disadvantages. It requires eight products to be 

10 formed, each of four values, and addition thereof. Hence it is not always best for |q 
high speed calculation. 

There is another disadvantage in this method. Although from one unit cube 
to the next the interpolated values are continuous, their derivative is not. That is, 
the slope of the interpolated values is discontinuous from one unit cube to the next, 

15 i.e. the line of the interpolated values bends sharply as we pass over the boundary. ^ 
Thus in practice a sharp step of color values will be apparent in the finished picture, 
and the cubic structure of the memory will show, to the detriment of quality. This 
effect can become quite serious. Fig. 4 shows a distribution of the interpolated values 
obtained according to the formula (I) which has a saddle form, which shows the 

20 aforementioned inconvenience clearly. An even continuous line of interpolated values 20 
in the unit square AjBidDx is obtained, and also in the unit square AjBjC^D^ but 
between these two squares, at their common border, the derivative of the interpolated 
values is discontinuous. . ■ 

The prior art interpolation methods, and their disadvantages, have been explamed 

25 above. A method according to the present invention will now be described. 25 
In Fig. 5, showing the two-rdimensional case, two adjacent interpolation regions 
AiBiCxDt and A 2 B 2 C 2 D 2 are shown. The centers of these unit squares are designated 
by Oi and O*, and interpolated values at these points are derived as averages of the 
function values at the four corners of the squares. Then interpolation is conducted 

30 linearly in each of the triangles AxC^B,, B1O1Q, dOiDj, DiOjAi, A 2 0 2 Bfi, B 2 0 9 Q, 20 
C 2 0 2 D 2 , and D 2 0 2 A 2 . That is, the point at which the value is to be interpolated is 
first checked to determine which of these triangles it falls into, and then the value at 
the point is determined by interpolation in the triangle in a fashion analogous to 
Fig. 2, by drawing lines from the point to the corners of the triangle, and then 

35 calculating the value of the function at the point as a weighted sura of the values 35 

at the corners of the triangle, giving each value at a corner a weighting of the ratio 
of the area of a second triangle whose corners are the point and the other two corners 
of the triangle, and the area of the triangle. In this method the magnitude of the 
discontinuity in the derivative of the interpolated values from one interpolation region 

40 to the next is much reduced. 40 
Now, considering the three-dimensional case, the basic interpolation method in 
a tetrahedral volume will be explained with respect to Fig. 6. Let ABCD be a tetra- 
hedron of which each vertex is a point at which the value of the function U to be 
interpolated is known. The value at point P, internal to the tetrahedron, is calculated 

45 as follows: draw lines from each vertex A, B, C, and D through the point P to meet 45 

the opposite sides of the tetrahedron in A', B', C, and D'. Then the interpolated 
value U(P) is 

U(A) X ratio of volumes of tetrahedra PBCD and ABCD 
+ 

U(B) X ratio of volumes of tetrahedra PDAC and ABCD 
+ 

50 U(C) X ratio of volumes of tetrahedra PDAB and ABCD 50 

+ 

U(D) X ratio of volumes of tetrahedra PABC and ABCD. 
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Now the ratio of the volumes of the tetrahedra PBCD and ABCD, for example, 
is the same as the ratio of the heights of P and of A from the plane of BCD. . 

Referring to Fig. 7, there is shown a unit cube interpolation volume ABCDEFGH, 
and the values of the function U of three variables are assumed to be known at the 

5 vertices of the cube. All three variations of the present method depend upon dissecting ^ 
this cube into tetraheldra whose vertices are either vertices of the cube, centers of 
faces of die cube, or the center of the cube. Then a series of comparisons are made 
to determine which of these tetrahedra contains the point at which the value of the 
function is required to be interpolated. Once this is determined, the value is then 

10 interpolated within that tetrahedron according to the method described above, using iq 
analytical geometry; It will be realised that .it is mathematically reasonable to 
interpolate, initially, the values of the function at centers of faces of the cube as the 
average of the values at the four corners of the faces, and the value of the function at 
the center of the cube as the average of the values at all eight vertices of the cube. 

15 Thus for each vertex of each tetrahedron of the dissection of Fig. 7 the value of the 15 
function is known, and therefore the method illustrated in Fig. 6 can be applied for 
interpolation/ 

In the dissection of Fig. 7 there exist twenty-four tetrahedra, of each of which 
one vertex is the center of the cubic unit, two vertices are vertices of the cubic unit 
20 which are connected by an edge of the cubic unit, and the fourth vertex is the center 20 
of a square face of the cubic unit one edge of which face is the said edge of the cubic 
. unit The twenty-four tetrahedra are all isomorphic. One of them is shown in Fig. 8. 

The centers of the faces of the cubic unit have been labelled as Qi Q 6 , and 

the center of the cube as O. The tetrahedron OABQ, is illustrated. The planes OAB, 
25 QiAB, OQA and OQ 2 A have the equations y f — z f =0, z f — 0, .x,+y f — 1=0, and 25 
. x,— y f '=0 respectively. Therefore it is dear that the condition for the point P to lie 
inside the tetrahedron OABQi is that y,-z f >0, x f +y f <0, and x f -yy>0. (Of 
course z f >0, by definition). Provided these conditions are all satisfied, the interpolated 
. value of the function may be calculated as outlined above. Thus U(P) is equal to 

U(A),X ratio of volumes of POQxB and OABQ, 

"4-' • ■ 

U(B) X ratio of volumes of POQxA and OABQi 

_|_ ' 

U(QO X ratio of volumes of POAB and OABQi 
•* 4* . 

U(O) X ratio of volumes of PQiAB arid OABQi 

=U(A) [1-*,-*] + U(B) [«-y f ] 

30 +U(Q0 [2(y,-z,] + U(O) [2 z,] .COT) ''. 30 

This is because the ratio of the volumes of the abovementioned tetrahedra, as - 
pointed out above, is the ratio of their heights, and the equations of their faces are 
. . . stated above. 

Similar results hold when the point P is in the other interpolation tetrahedra. 
35 A complete table of the conditions for discrimination of which tetrahedron contains 35 

the point P, and of the factors which are used for calculation of the interpolated value 
in each case, is shown as Table 1. Using this table, by testing the conditions that 
are not parenthesized, it is possible to characterise the tetrahedron which contains 
the point P, and accordingly it is not necessary to test the parenthesized conditions. 
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CALCULATING FACTORS 
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It is readily understood that the calculation is far simpler in practice than the 
method of the abovementioned formula (II). Further, in this method, the dis- 
continuities across the borders between one unit cube and the next are much reduced, 
since the values near the face of the unit cube are much more dominated by the 
5 values at the four comers of the face than in the prior art method of (II). 5 
In fact die color picture output signals in the memory commonly vary mono- 
tonically, and therefore a more simple and coarse method of interpolation than the 
one outlined above may well be satisfactory in a particular case. Therefore the method 
of Figs 9 and 10 may well be acceptable, although it is not quite so accurate as the 
10 method of formula (III). In Fig. 9 is shown a dissection of the unit cube into 10 
tetrahedra all of those vertices are vertices of die unit cube. Thus this method nas 
the advantage that no averaging of values at the vertices of the unit cube is necessary 
in order to determine values at the centers of the faces of the unit cube and at its 
center. 

15 The unit cube is dissected into six tetrahedra by three planes which have , a line 15 

J in common which is the long diagonal of the unit cube, and each plane is inclined 
to the other two at 60° and contains two edges and four vertices of the unit cube. A 
typical one of the six tetrahedra is illustrated in Fig. 10. In this case the conditions 
for the point P to lie within this tetrahedron are that x,>y,>z„ as can be easily . 

20 worked out using solid geometry, as before. In the same way the interpolated value 20 
U(P) is equal to 

U(A) [l-x f ) + U(B) [x,-y f ] +U(C) [y t -Zt] + U(D) z,. 

Similar discriminating conditions and calculating factors can be worked out for 
the other five tetrahedra. Table 2 shows the complete set It is readily appreciated 
25 that this variation of the method is easier in calculation than the method of formula 25 
(HI), albeit at a slight loss in accuracy. 
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of the method, illustrated in Fig. 7. However, this variation will only be slight in the 
case of monotonic functions, and is quite tolerable. 

WHAT WE CLAIM IS:— 

1. A linear interpolating method for color signals in a memory of a picture 
5 reproducing machine, comprising storing appropriate values of color picture output 

signals corresponding to certain stepped values of color input signals in the memory 
addressed in a three-dimensional fashion, and interpolating values of color output 
signals at points which are between said values by: 

dividing up the cubic interpolation unit of the memory which is constituted by 
10 a single step of each of the color input signals into a plurality of tetrahedra whose 
vertices are either vertices of the cubic unit, centers of its faces, or its center point; 

calculating the color output signal at each vertex of these tetrahedra which is a 
center of a face of the cubic unit, if any, by averaging the values of the color output 
signal at the four vertices which are corners of said face, and at the center point of 
15 the cubic unit by averaging the values of the color output signal at all eight of the 15 
vertices of the cubic unit; 

determining which of these tetrahedra includes the interpolation point at which 
the value of the color output signal is to be interpolated ; and 

deriving the interpolated value at the interpolation point as a weighted sum of 
20 the values at the four vertices of the determined tetrahedron, the value at each vertex 20 
being given a weight corresponding to the ratio of the volume of a second tetrahedron 
whose vertices are the interpolation point and the other three vertices of the deter- 
mined tetrahedron to that of the determined tetrahedron. 

2. A method as in Claim 1, wherein the cubic unit is divided into twenty-four 

25 tetrahedra, of each of which one vertex is the center of the cubic unit, two vertices 25 
arc vertices of the cubic unit which are connected by an edge of the cubic unit, and 
the fourth vertex is the center of a square face of the cubic unit one edge of which 
face is the said edge of the cubic unit. 

3. A method as in Claim 1, wherein the cubic unit is divided into six tetrahedra 

3® by three planes which have a line in common, said line being a long diagonal of the 30 
cubic unit, and each plane containing two edges and four vertices of the cubic unit 

4. A method as in Claim 1, wherein the cubic unit is divided into five tetrahedra 
by four planes, each of which contains exactly three vertices of the cubic unit; said 
planes intersecting one another along lines which are diagonals of faces of the cubic 

35 unit 35 

5. A linear interpolating method for colour signals in a memory of a picture 
reproducing machine as claimed in Claim 1 and substantially as hereinbefore described 
with reference to and as illustrated in Figs. 5, 6, 7 and 8, 9 and 10 or 11 and 12 of. 
the accompanying drawings. 

H. L, COTTRELL & CO., 
Chartered Patent Agents, 
22 Whitefriargate, 
Hull, HU1 2EX. 
Agents for the Applicants. 
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